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There are indications that hadronic loops in some electroweak observables are almost saturated 
by parton level effects. Taking this as the hypothesis for this work, we propose a genuine parton 
level estimate of the hadronic light- by- light contribution to the anomalous magnetic moment of the 
muon, a^ BL (had). Our quark mass definitions and values are motivated in detail, and the simplicity 
of our approach allows for a transparent error estimate. For infinitely heavy quarks our treatment is 
exact, while for asymptotically small quark masses aJ; BL (had) is overestimated. Interpolating, this 
suggests to quote an upper bound. We obtain aJ; BL (had) < 1.59 x 10 -9 (95% CL). 

PACS numbers: 13.40.Em, 14.60.Ef, 12.20.Ds, 12.38.Bx. 



The E-821 Collaboration at BNL has measured 
the anomalous magnetic moment of the muon, a M = 
(gy, — 2)/2, with an uncertainty of Aa^ = ±0.63 x 10~ 9 . 
This would provide sensitivity to new physics scales, A, 
as large as A ~ mu/^Asj, = 4.2 TeV, where is the 
mass of the muon 0. A new experiment at BNL is ap- 
proved aiming at even greater (±0.2 ppm) precision and 
implying a reach up to 7 TeV. Unfortunately, the inter- 
pretation of a M is compromised by large theoretical uncer- 
tainties introduced by hadronic effects diluting the new 
physics sensitivity. These mainly arise from two-loop vac- 
uum polarization (VP) effects, a^ p (2, had), and from the 
three-loop contribution of light-by-light type, aJ; BL (had). 
The calculations of a^ BL (had) based on chiral perturba- 
tion theory are the only ones to date solidly based 
on a systematic expansion. However, the estimated un- 
certainty, > ±10~ 9 , is significantly larger Q than the 
measurement error. There is a variety of model esti- 
mates, all supplementing the dominant 7r°-exchange con- 
tribution with other resonance exchanges and ^-loops. 
Current analyzes 00 agree reasonably well on the mag- 
nitude (where residual differences are largely understood) 
and the sign of a^ BL (had) . They have reached a high level 
of sophistication, but the error estimates remain rough 
guesses. There is also an estimate based on the in- 
stanton liquid model with a small (< 10%) quoted error. 

In this work we offer a complementary way to estimate 
a^ BL (had), with no need to commit to the dominance of 
any particular type of contribution. It is a naive par- 
ton level estimate which is solid in the heavy quark limit 
where it matches onto perturbative QCD, but overesti- 
mates the light-by-light contribution in the chiral limit. 
Therefore, our approach naturally implies an upper bound 
for a^ BL (had), which is very useful in view of the a M mea- 
surement lying above the Standard Model prediction. It 
can also be applied to a^ p (had) , which serves as a refer- 
ence case and allows for a transparent error estimate. 

Attempting to obtain aJ; BL (had) directly at the par- 
ton level seems hopeless at first since perturbative QCD 
cannot be applied except for the heavy c and b quarks. 



Moreover, in the chiral limit, in which the bulk of the 
effect arises from light, charged, and essentially point- 
like pseudo-Goldstone bosons, one has a^ BL (had) < 
from scalar QED 0, while leptons and heavy quarks 
contribute positively However, as shown in Ref. 0, 
the 7r ± -loops receive large chiral symmetry breaking cor- 
rections, rendering the net -^-contribution almost negli- 
gible. Thus, at least as far as aJ; BL (had) is concerned, 
we are dealing with a kinematic regime that conceiv- 
ably admits an alternative description in terms of quark 
in place of chiral degrees of freedom. Indeed, an esti- 
mate of the typical momenta, (j,, in the 7r ± -loops yields 
/i ~ 4.25 w 0.6 GeV, i.e. large enough to feel indi- 
vidual partons even within pions. These observations 
revive the question whether hadronic light-by-light dia- 
grams are best described by hadrons or by partons or 
a mix. We are aware of three parton based attempts 
in the literature. An early approach using exclu- 
sively constituent quarks [12j resulted in a numerical es- 
timate, albeit with the wrong sign. Ref. [lj| opted for 
a mixed treatment in which constituent quark masses, 
m u d = 0.33 GeV, m s = 0.5 GeV, and m c = 1.87 GeV 
(with assumed errors of ±10%), were employed in the an- 
alytical result of Ref. |l(J and produced the correct sign. 
The adopted definition, values, and uncertainty were not 
further justified, but in any case the constituent quarks 
were assumed to yield only the high momentum contri- 
bution which was supplemented with 7r°-exchange graphs 
as described by the chiral constituent quark model 14] . 
Ref. 0] is very close to both our strategy and philosophy 
and will be commented on below relation (|10|> . 

Indeed, a main obstacle for any model involving quark 
degrees of freedom at hadronic scales is to find the appro- 
priate quark mass definition. In contrast to pole masses, 
MS masses, m q (^i) (where fj, is the renormalization scale), 
tend to produce small perturbation coefficients (at least 
for sufficiently inclusive quantities). We will therefore 
avoid pole masses for heavy quarks and similarly con- 
stituent masses for light quarks. The scale choice is 
implied by the requirement of avoiding spurious loga- 
rithms, ln/^ 2 /m^(/i) or ln^ 2 /m 2 . They can be nullified 
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(resummed) by choosing fi = m q (fJ,) for m q (rh q ) > 
and n = for m q (rh q ) < m^, i.e. the higher of the two 
masses sets the scale. There are also genuine logarithms 
of the form In m q (/x) / mz which cannot be eliminated. 

But for the light (u, d, s) quark contribution to — 2, 
the values of m q {rh q ) or rh q (m^) are unknown because 
d s and consequently the anomalous quark dimensions 
formally diverge. In a similar context, however, it was 
possible to proceed phenomenologically. Namely, the 
hadronic contribution to the renormalization group evo- 
lution (RGE) of the QED coupling (the running), a(n), 
can be mimicked |l6j by effective quark masses. These 
have been defined more precisely in Ref. |17| as threshold 
quark masses, fh q , and correspond to the scale at which 
the theory is changed to include or exclude the quark 
q and where the one-loop QCD (3 function coefficient 
changes correspondingly. Since this is not a perturba- 
tive description, higher order QCD effects are assumed 
to be absorbed in the threshold masses, leading to |l7|. 



rh u = fh d = 176 ± 9 MeV, 

fh s = 305^ MeV, 

m c = 1.18 ± 0.04 GeV, 

fh b = 4.00 ± 0.03 GeV, 



(1) 
(2) 
(3) 
(4) 



where isospin symmetry, m u {m u ) — rhd{rhd), was im- 
posed. The uncertainties in the first (last) two lines have 
a correlation of -100% (+29%). The values (JTHUl are 
implied by the considerations [17J in the next paragraph. 

Heavy quark masses are qualitatively given by one 
half of the mass, Mis, of the corresponding IS quarko- 
nium resonance. If the quark is heavy enough that 
QCD perturbation theory can be applied, this state- 
ment can be made quantitative resulting in well stud- 
ied QCD sum rules which are dominated by the 
lS'-quarkonium resonance and supplemented by higher 
vector resonances and a continuum contribution. They 
yield expressions for rh q (m q ) with small QCD corrections 
which are known |2flj | and included up to order a 2 .. The 
threshold masses in Eqs. iJSHlJl can then be obtained in 
terms of these using Eq. (32) in Ref. 01- Defining the 
parameter, £ 9 = 2fh q jM\s |l7|, we have that asymptot- 
ically £ g — > 1 for fh q — > co and £ g — > for fh q — > 0. 
Thus, for a heavy quark £ g ~ 1, while for a light quark 

« 1. Also, we expect fh\ < ffi2 => £i < £2, which 
can be checked explicitly by comparing = 0.845±0.006 
with £ c = 0.76 ± 0.03. This implies the upper bound 
£, s < £,c- A lower bound, £ s > 0.47, can be obtained by 
considering the SU(3) limit, £ u = £ d = £ s . Since the 
total contribution to the runnin g of a due to u, d, and s 
quarks is known experimentally 34] , this (together with 
isospin symmetry) implies Eq. J5J ■ Assuming in addition 
the supp ression of flavor singlet contributions (the OZI 
rule |3Jj) then implies Eq. JJJ and explains the complete 
anti-correlation with Eq. (j2J). Isospin violations will, of 
course, induce a small mass difference, fhd — fh u . 



In Ref. [13 the threshold masses were used for 

the RGE of the weak mixing angle, sin 2 #(/i), which is 
rigorously justified as the same vector current correlator 
convoluted by the same weight function enters in both, 
d(/x) and sin 2 #(/i). What we are proposing here, is that 
a^ BL (had) can be modeled by treating quarks like massive 
leptons and using the values JBHJl. Notice, that these 
values are systematically lower than typical constituent 
masses |l2|, which is reasonable since we do not add an 
independent contribution from pion degrees of freedom. 

It is a priori not obvious that this procedure and 
these values for the quark masses (in fact, any set of 
values) will generate the correct a^ BL (had), but this rep- 
resents the only ad hoc model assumption that we need 
to make. One justification is the observation mentioned 
before, that the relevant kinematic region conceivably ad- 
mits a description in terms of quark degrees of freedom, 
which is reinforced by the lower bounds, £ s > 0.47 and 
£ Ut d > 0.42, showing that the physical light quark masses 
are midway between the perturbative (£ — > 1) and chi- 
ral (£ — > 0) limits. Thus, at a qualitative level, both 
a^ BL (had) and the RGE of d(/x) appear to be describ- 
able by moderately heavy quark degrees of freedom. 

Although our treatment clearly improves with growing 
fhq and becomes exact for fh q — > co, it is important to 
establish quantitative tests and to apply our approach to 
physical quantities for which the answer is known with 
sufficient precision. For this we note that a^ p (2,had) is 
described by the same vector current correlator as the 
RGE of a(fi) but is convoluted with a different weight 
function, K (1 ^(s) 21], so that just as in the case of 
a^ BL (had) it is a priori not obvious that our approach 
and quark mass values will generate the correct answer, 

a^ p (2,had) x 10 9 = 70.30 ± 0.43 ± 0.88 = 70.3 ± 1.0 (5) 

To arrive at Eq. J^J we averaged the results based on 
e + e~ annihilation data and r decay data which are both 
taken from Ref. 22] . The second uncertainty accounts for 
the discrepancy between the two data sets. It is chosen 
to one half of the difference in central values and added in 
quadrature to the first (experimental) one. We compare 
this result with the analytical expression [23I l24j , 



a vP( 2 ,had) = ^iVc^Q^ 

q n— 



2n+2 



(6) 



nln(m 2 /m 2 ) 



8n A + 28n 2 - 45 



(n + 3)(2n + 3)(2n + 5) [(n + 3)(2n + 3)(2n + 5)p 



yielding, 



a^ p (2, had) = 65.1+^ T 0.1 x 10~ 9 , 



(7) 



The first and second uncertainty are induced, respec- 
tively, from Eqs. and Eqs. (OQ. Effects due to 
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raid + 1 m u cancel to first order in a^ p (2,had) (but not 



-,LBL / 



,, "(had) below) and can be neglected. Notice, that 

the estimate is in reasonable agreement with Eq. (J5J 
and that our central value is below it. This is consis- 
tent with the chiral limit, — > with m^/m^ fixed, 



underestimating 



(2, had). Indeed, in this limit, 



0.39 



(where \x is a fixed reference scale) so that m^/fhq — > 0. 
Combined with K(s) > this implies a systematic un- 
derestimate. The fact that the estimate Q reproduces 
the experimental result within about 10% provides phe- 
nomenological support for our approach and its use of 
quark degrees of freedom. We can also directly compare 
the contribution from c quarks, a^ p (2,c) = 1.27 x 10~ 9 
from Eqs. and 0, with the analytical result up to 
order a 2 s originally obtained in the second Ref. [20| in 
the form given in Eq. (4) of Ref. [25J. Using m c (rh c ) — 
1.29 GeV and a s (M z ) = 0.1216 ± 0.0017 (correspond- 
ing to a s (m c ) = 0.432 and a s (fh c ) — 0.473), the latter 
yields a^ p (2,c) = 1.46 x 10~ 9 . Thus, here we underes- 
timate the known result by 13%, which is traceable to 
different a s dependences in the RGE of a(/j,) and the 
leading 0(mf 1 /ml) term in a^ p (2,had). QCD correc- 
tions to the C(m^/m^) term which are significant only 
for lighter quarks tend to have a compensating effect, of- 
fering a perturbative rationale why our approach seems 
to work better for light quarks than for charm. 

The success of the test in the previous paragraph may 
be coincidental. Therefore, we performed three similar 
tests on various three-loop VP effects, aY p (3i, had), each 



additional electron (hadron) loop insertion. Our method 
(in parentheses) reproduces the central values of Ref. [2T 



,VP/ 



'(3a, had) = -2.11 x 10" 9 (-1.82 x 10" 9 ), 
5 (3b,had)= +1.07 x 10~ 9 (+0.99 x 10~ 9 ), 
a^ p (3c,had)= +2.7 x 10~ n (+2.8 x lO" 11 ), 



within 14%, 7%, and 4%, respectively. Note, that these 
represent genuine "hit or miss" tests involving no ad- 
justable parameters, and all four passed at the 10% level. 

We finally apply our approach to alt BL (had). Our mas- 
ter formula is the one for heavy leptons 10], 



LBL 



3 00 2n 2 



my, 

n=l 1 k=0 



where Ci,o = 3/2( 3 — 19/16, Ci : i = — 0, and where 
the other C n ,k up to n = 5 can be found in Eq. (5) of 
Ref. 0. The central values in Eqs. then imply 



is the model error, determined by the 9.5% difference 
between the upper end of the range (0 from the cen- 
tral value in Eq. JTJl. Similarly, the upper and lower er- 



ror in Eqs. and 10 correspond to 



.LBL 



(had) = 



1.27 ± 0.18 x 10~ 9 and aJ; BL (had) = 1.47 ± 0.04 x 10~ 9 , 
i.e. model errors of 14.1% and 2.9%, respectively. As 
expected, b quark effects are negligible, but the contri- 
bution from c quarks (usually ignored) turns out to be 
0.04 x 10 -9 . Although small it might serve as an interest- 
ing testing tool of other model calculations of a^ BL (had) 
which should be consistent with what we have obtained, 
since our method is solid (up to effects of 0(a s )) for 
heavy quarks. To account for isospin breaking we shift 
fn-d by +2 MeV and rh u by —1 MeV, increasing all results 
by 0.01 x 10~ 9 to 



LBL 



(had) = 1.37; 



-0.27 
-0.15 



X 10" 



(9) 



with a different kernel function, K^l-s) 26, 27] . i = a 
contains an additional photonic correction or muon-loop 
relative to ajf p (2,had), while i = b (c) corresponds to an a LBL (had) 



where the error covers the three ranges above and is 
constructed to include the intrinsic model uncertainty. 
As in any hadronic model there would be additional 
model uncertainty if the comparison to other cases (here 
a^ p (had)) was regarded as pure coincidence. The re- 
sult, 1.36 ± 0.25 x 10~ 9 , of Ref. is the first to take 
short-distance QCD constraints explicitly into account. 
By construction, this is a feature of our approach, as 
well, so it is gratifying that the estimate © turns out to 
agree perfectly. Our central value is higher than previ- 
ous ones 0, which is expected given that in the chiral 
limit, — > with m^/m^ fixed (and presumably also 
with m M fixed, but the scalar QED result in this case is 
not known), we overestimate a^ BL (had) since rr^-loops 
contribute negatively. 

Now we repeat the above analysis with the model 
errors multiplied by 1.645 to shift from uncertainties 
estimated in a "ler spirit" to 90% ranges, and find 



1.37+0.21 xl0~ 9 , aJ; BL (had) 



1.28±0.30x 



10~ 9 , and aJ; BL (had) = 1.48 ± 0.07 x 10~ 9 , respectively. 
The upper errors have converged, and since as mentioned 
our approach tends to overestimate a^ BL (had), we quote 
as our final result the 95% CL upper bound, 



LBL 



(had) < 1.59 x 10" 



(10) 



,LBL 



In a similar paper ,15], a 1-parameter fit to (2, had) 
resulted in m u = = m s — 180 MeV = 166 ± 1 MeV. 
These values are within errors consistent with Eqs. 00, 
reflecting our statement that a^ p (2,had) and a(fj,) can 
be approximated using the same masses. We benefited 
from a deeper m^/m^-expansion of a ^ BL Q, lowering the 
light quark contribution relative to |l5j by ~ 6%. The 
most important new aspect here is that the recent scalar 
QED calculation allows us to interpolate between the 
chiral and heavy quark limits. If one can prove (conceiv- 
ably by lattice simulations) the monotony of the product 
(had) as a function of m q , our bound (|10J) would 



-,2_,LBL 



(had) = 1.36 ± 0.13 x 10 , where the uncertainty become rigorous, but this remains a loophole for now. 
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We stress that it is worthwhile to compute the a s cor- 
rection to the heavy quark light-by-light contribution to 
to check explicitly its magnitude and more impor- 
tantly its sign, to compare it to the RGE of a(n), and 
to use it as a refinement of our method. This contribu- 
tion may be called A^ 2 ^ {m^/rh q ) and corresponds to a 
subclass of 60 of the 469 Feynman diagrams contribut- 
ing at 0(a ) to what is known as A 2 (m^/m T ) with m T 
replaced by rh q (m q ) and one prefactor of a by a s (rh„). 
This is possible with state-of-the-art computer codes [29j . 

We have developed (quark mass definition, scale set- 
tings, quark mass values, and uncertainties) what we be- 
lieve are necessary ingredients for a parton model esti- 
mate of aJ; BL (had) and similar quantities. The model has 
four key benefits: (i) it is simple and transparent, (ii) it 
can be tested, as we did successfully using a^ p (had), (iii) 
short-distance constraints are taken into account by con- 
struction, and (iv) it naturally provides an upper bound 
on a^ BL (had) because, as we have argued, our method 
tends to yield an overestimate and also because the 95% 
upper error is virtually insensitive to our input masses. 

We would like to compare our approach with cele- 
brated models of low-e nerg y QCD, such as the vector me- 
son dominance model |30| . the \jNc expansion [3l|], the 
chiral constituent quark model 01 1 or the extended 
Nambu-Jona-Lasinio [13 model. It is not obvious a 
priori that approximations like Nc = 3 > 1 |3ll ] or 
Aqcd ^ My; b (the spontaneous chiral symmetry break- 
ing scale) [l4| . or the dominance of certain resonances [3(J 
or operators — while plausible — would generate 
useful results. But experience showed that they do. The 
assumption made in this work is that for the type of ob- 
servables considered, parton level effects dominate and 
almost saturate hadronic loops. It is in line with the 
observations of Ref. Q, and we found further evidence 
for it in both, a^ p (had) and a^ BL (had), which passed 
quantitative tests with better than expected phenomeno- 
logical success. This may be interpreted as pure coinci- 
dence and, clearly, it is important to submit our method 
to other case studies in the future. More optimistically, 
however, our results offer an alternative view on hadronic 
loops and may lead to a rigorous bound on a^ BL (had). 
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